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We present a calculation of the N to A electromagnetic transition amplitudes using the method 
of QCD sum rules. A complete set of QCD sum rules is derived for the entire family of transitions 
from the baryon octet to decuplet. They are analyzed in conjunction with the corresponding mass 
sum rules using a Monte-Carlo-based analysis procedure. The performance of each of the sum rules 
is examined using the criteria of OPE convergence and ground-state dominance, along with the 
role of the transitions in intermediate states. Individual contributions from the u, d and s quarks 
are isolated and their implications in the underlying dynamics are explored. Valid sum rules are 
identified and their predictions are obtained. The results are compared with experiment and other 
calculations. 
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I. INTRODUCTION 

The determination of low-energy electromagnetic prop- 
erties of baryons, such as charge radii, magnetic and 
quadrupole moments, has long been a subject of interest 
in nucleon structure studies. For example, the E2/M1 ra- 
tio of the transition amplitudes in the process "/N — > A is 
a signature of the deviation of the nucleon from spherical 
symmetry. 

The first serious attempt to determine the ratio from 
quantum chromodynamics (QCD), the fundamental the- 
ory of the strong interaction, was the lattice calculation 
of Leinweber et al. [l| . The calculation has been updated 
recently by Alexandrou et al. 0, Q with better technol- 
ogy and statistics. In this work, we present a calcula- 
tion of the transition amplitudes using the SVZ method 
of QCD sum rules [3] which is another nonpcrturbative 
approach firmly-entrenched in QCD with minimal mod- 
eling. The approach provides a general way of linking 
hadron phenomenology with the interactions of quarks 
and gluons via only a few parameters: the QCD vacuum 
condensates and susceptibilities. The analytic nature of 
the approach offers an unique perspective on how the 
properties of hadrons arise from nonperturbative inter- 
actions in the QCD vacuum, and gives a complimentary 
view of the same physics to the numerical approach of 
lattice QCD. It has been successfully applied to almost 
every aspect of strong-interaction physics, including the 
only other attempt on the N to A transition that we 
know of by lofFe [5| over 25 years ago. It was a limited 
study conducted in the context of the magnetic moments 
of the nucleon, and the results were inconclusive due to 
large contamination from the excited-state contributions. 

The above discussion only concerns with the classic 
QCD sum rule method by SVZ. There are a number of 
studies on N to A transitions in the so-called light cone 
QCD sum rule method (or LCSR) [6, 7]. The main dif- 
ference is that it employs as input light-cone wavefunc- 
tions, instead of vacuum condensates. One advantage of 
the light-cone QCD method is that one can compute the 
transition form factor at non-zero momentum transfer 
relatively easily [1, [§] . 



Our goal is to carry out a comprehensive study of the 
N to A transition in the SVZ QCD sum rule method with 
a number of features. First, we employ generalized inter- 
polating fields which allow us to use the optimal mixing 
of interpolating fields to achieve the best match in the 
QCD sum rules. The interpolating field used previously 
is a special case. Second, we derive the complete tensor 
structure and construct a new set of QCD sum rules. The 
previous study was quite limited in its tensor structure. 
Third, we study all tranisitons from the octet to te deu- 
plet, not just the proton channel. Fourth, we perform a 
Monte-Carlo uncertainty analysis which has become the 
standard for evaluating errors in the QCD sum rule ap- 
proach. The advantage of such an analysis is explained 
later. Fifth, we use a different procedure to extract the 
transition amplitudes and to treat the transition terms 
in the intermediate states. Our results show that these 
transitions cannot be simply ignored. Sixth, we isolate 
the individual quark contributions to the transition am- 
plitudes and discuss their implications in the underlying 
quark-gluon dynamics. A similar study has been per- 
formed for the magnetic moments of octet baryons [lO|. 



The paper is organized as follows. In Section |TT1 the 
physics of N to A transition is reviewed. In Section IIIIl 
the external field method of QCD sum rules is applied 
to this process. Master formulas are calculated using the 
full sets of interpolating fields for the octet and decuplet 
baryons. Then the full phenomenological representation 
and QCD side are calculated, leading to a complete set 
of QCD sum rules for the entire family of transitions 
from the octet to the decuplet. Section HVl describes our 
numerical procedure to extract the transition amplitudes. 
Section |V] summarizes the results and gives a comparison 
of our results with experiment and other calculations. 
The conclusions are given in Section IVII 
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II. PHYSICS OF N TO A TRANSITIONS 
A. Definition of the form factors 

The concept of form factors plays an important role 
in the study of internal structure of composite particles. 
The non-trivial dependence of form factors on the mo- 
mentum transfer {i.e., its deviation from the constant 
behavior) is a signal of the non-elementary nature of the 
investigated particle. The current matrix element for the 
jN — !■ A transition can be written as 

<A{k',s')\j^\N{k,s) >^uf{k',s')rp^u{k,s), (1) 

where k, s and fc', s' denote initial and final momenta 
and spins, the electromagnetic current, u the spin- 
1/2 nucleon spinor, and the Rarita-Schwinger spinor 
for the spin-3/2 A. This matrix element is the most 
general form required for describing on-shell nucleon and 
A states with both real and virtual photons. There exist 
different definitions in the literature for the form factors 
describing the vertex function Fp^. One is in terms of 



Gi, G2, and G3 by Jones-Scadron [11|: 



where the invariant structures are 



PA* 



9 ■ 

- q 



75pp; 
-Pffpp; 



9^5pp- 



(2) 



(3) 



Here q — k' — k is the four-momentum transfer. For real 
photons (g^ = 0), only Gi(0) and G2(0) play a role. Note 
that the form factors defined this way are dimensionful: 
Gi in GeV-\ G2 in GeV-^, and G3 in GeV^^ xhey are 
analogous to the Pauli-Dirac form factors Fi and F2 for 
the nucleon. 

Another definition commonly used in lattice QCD cal- 
culations [il, 0] is 



<A{k',s')\j'^\Nik,^s)> 



(4) 



where 



F"^ = GMl{q')K^, + GE2{q^)KZ + 



Gc2{q^Wc2, 



(5) 



in terms of the magnetic dipolc Gmi, electric quadrupole 
Ge2, and Coulomb quadrupole Gc2 form factors. They 
are analogous to the Sachs form factors Ge and Gm for 
the nucleon. The invariant functions involve the masses 
of the particles explicitly, 



"■Ml 
r^PP 
"E2 



"C2 
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where P = [k' + k)/2 and Vl[q^) = [(toa + tunY - 

We stress that the two sets of form factors are equiva- 
lent in describing the physics of N to A transitions. They 
are related by 
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Here we have switched to the notation of = —q^ since 
most studies of form factors are for space- like momentum 
transfers (g^ < or > 0). The inverse form, omitting 
the explicit dependence in the G's, is given as: 



Gi = 
G2 = 

G3 = 



-3(Ge2 - Gmi)[toa("ia + ruN)] 
/[2mN{{mA + rriNf + Q^)]; 
-3{mA + mN)[i-2Gc2^^ 
+GMi{m\ - 2mAmN - 
+G_E2(— 3m^ -I- 2mAmN -r m-j^ 
/[2mN{(m\ - TO|r - (5^)2 ^ 4^ 



m 



N 
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-3m| 

*A ~ "'TV ^ ^ ) ^'*"'-Ar 

-3(toa + TOAr)[(-2Gc2(m-A ^ "^w) 
-|-Gmi(wa - 2mAmN + rn%, + Q^) 
+GE2{^rn\ + '^rUAmN + m% + Q'^))] 
/[4m^((mi - m% - Q^)^ + im%Q^)] 



(8) 



From these relations, we see that G^f 1 is on the same or- 
der as Gi, Ge2 is proportional to Gi/mA + G2 at = 0. 
The smallness of G£;2 comes from cancelation between Gi 
and G2, making it quite difficult to quantify accurately. 

The Particle Data Group (PDG) [ij uses the transi- 
tion amplitudes /mi and fE2 which are re-scaled versions 
of the Sachs form factors 



/a/1 

fE2 
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-Ge2, 



where e — V Ana. In the rest frame of the A at g^ = 0, 
energy-momentum conservation sets 2|(7|mA = m'^ — 1 
so the relations simplify to 



'AT 



J Ml = 2;;r77'- 2m« ) Gmi, 

fp2 = ^/mi-m^.i/2g 

J^^ 2mN^ 2mM > ^E2- 



(10) 



They are related to the well-known helicity amplitudes 

by 



= 5^(3^3/2 + 73^1/2); 

fE2 = 573(^3/2 - %/3^1/2)- 

A commonly quoted quantity is the ratio 

/b2 Ge2 



Rem 
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Gmi 
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Partial photonuclear decay widths may also be related 
to the Sachs form factors assuming continuum dispersion 
relations, 



^E2 = 



3q (mj-mw) 



^2 

"-"Ml' 



(13) 



B. Experimental information 

Throughout this work, we use the generic term N to 
A transition to refer to the eight transitions from the 
baryon octet to the decuplet: 



A+, 
AO, 

E*- 

"*0 



(14) 



IP 
7n 
7S+ 

7S- 
7SO 
7^ 
7A 



Experimentally, pion photo- and or electro-production 
has been used to access the transition amplitudes Gm 
and Ge- The only transition measured so far is 7p — > 
A+. The Laser Electron Gamma Source (LEGS) at 
Brookhaven has made measurement of proton cross sec- 
tions and photon asymmetries. Experiments on pion 
electroproduction were performed at the Mainz Mi- 
crotron (MAMI) in the A resonance region and low mo- 
mentum transfer. A program using the CEBAF Large 
Acceptance Spectrometer (CLAS) at Jefferson Lab has 
been inaugurated to improve the systematic and statisti- 
cal precision by covering a wide kinematic range of four- 
momentum transfer Q^. At the MIT-Bates linear accel- 
erator, an extensive program has also been developed. 



TABLE I: A summary of experimental results for the ratio 
Rem in the 7p transition. 



Experiment 




Gmi 


Ge2 


Rem 




(GeV") 




{^J.N) 


(%) 


MAMI [13] 


0.06 


2.51 


0.057(14) 


-2.28(62) 


MAMI [14] 


0.38 


3.01(1) 


0.096(9) 


-3.2(0.25) 


LEGS fl5] 









-2.5 


JLab fl6] 


0.4 


1.16 


0.040(10) 


-3.4(4)(4) 


Bates [17] 


0.126 


2.2 


0.050(9) 


-2.0(2)(2) 


PDG fl2] 





3.02 


0.050 


-2.5(0.5) 



Table U lists some of the main experimental data and 
analysis results at low momentum transfer. It is clear 
from the results on Rem that the quadruple component 
in the N to A transition is quite small, signaling a slight 
deformation of the nucleon. The minus sign implies that 
its shape is oblate rather than prolate. The bulk of the 



transition is dominated by the magnetic dipole compo- 
nent, which results from a spin-flip in one of the quarks. 
The PDG results corresponding to Gmi and Ge2 in other 
notations are Gi = 2.70 GeY-\ G2 = -1.65 GeV"^; 



/2 



0.135 GeV-1/2, A3/2 = -0.250 GeV-1/2. 



-1/2. 



/mi ^ 0.284 GeV^^/^ fE2 = -0.0047 GeV" 

The extraction of the small Ge2 amplitude is very 
sensitive to the details of the models and the specific 
database used in the fit. Three commonly used are the 
dynamical model of Sato and Lee (SL model) fl8|, the 
unitary isobar model MAID [13], and the partial- wave 
analysis model SAID [2(| . The SL model uses an effective 
Hamiltonian consisting of bare vertex interactions and 
energy-independent meson-exchange ^ A transition 
operators which are derived by applying a unitary trans- 
formation to a model Lagrangian describing interactions 
among 7, 11, p, uj, N and A fields. With appropriate phe- 
nomenological form factors and coupling constants for p 
and A, the model gives a good description of ttN scatter- 
ing phase shifts up to the A excitation energy region. The 
MAID and SAID models are partial-wave-based models 
for pion photo and electroproduction. It start from the 
Lagrangian for the nonresonant terms, including explicit 
nucleon and light meson degrees of freedom coupled to 
the electromagnetic field. The resonance contributions 
are included by taking into account unitarity to provide 
the correct phases of the pion photoproduction multi- 
poles. This model provides a good description for indi- 
vidual multipoles as well as differential cross sections and 
polarization observables. 



III. QCD SUM RULE METHOD 

The starting point is the time-ordered correlation func- 
tion between a nucleon state and a A state in the QCD 
vacuum in the presence of a static background electro- 
magnetic field Ffj^^: 

n«(p) = z Jd^xe^P^O\T{f^^{x)fj''{0)}\Q)E, (15) 

where r]^ {a Lorentz index) and ry^ are interpolating field 
operators carrying the quantum numbers of the baryons 
under consideration. The task is to evaluate this correla- 
tion on two different levels. On the quark level, the corre- 
lation function describes a hadron as quarks and gluons 
interacting in the QCD vacuum. On the phenomeno- 
logical level, it is saturated by a tower of hadronic in- 
termediate states with the same quantum numbers. By 
matching the two, a link can be established between a 
description in terms of hadronic degrees of freedom and 
one based on the underlying quark and gluon degrees of 
freedom as governed by QCD. The subscript F means 
that the correlation function is to be evaluated with an 
electromagnetic interaction term. 



(16) 
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added to the QCD Lagrangian. Here ^4^ is the external 
electromagnetic potential and = Cqq^f^q is the quark 
electromagnetic current. 

Since the external field can be made arbitrarily small, 
one can expand the correlation function 

n„(p) = nW(p) + nW(p) + ---, (i7) 

where 11^°^ (p) is the correlation function in the absence 
of the field, and gives rise to the mass sum rules of the 
hadron. The transition amplitudes will be extracted from 
the linear response function !!„ ■* (p) . The action of the ex- 
ternal electromagnetic field is two-fold: on the one hand 
it couples directly to the quarks in the hadron interpo- 
lating fields; on the other it polarizes the QCD vacuum. 
The latter can be described by introducing new parame- 
ters called QCD vacuum susceptibilities. 



Here u, d, and s are up, down, and strange quark field 
operators, C is the charge conjugation operator, the su- 
perscript T means transpose, and £abc ensures that the 
constructed baryons are color-singlet. The real parame- 
ter P allows for the mixture of the two independent cur- 
rents. The choice advocated by loffe [21| and often used 
in QCD sum rules studies corresponds to /3 = — 1. In 
this work, /3 is allowed to vary in order to achieve max- 
imal overlap with the state in question for a particular 
sum rule. The normalization factors are chosen so that in 
the limit of SU(3)-flavor symmetry all N to N correlation 
functions simplify to that of the proton. 
For the spin-3/2 baryon decuplet, we use 



A. Interpolating fields 



We need interpolating fields for the baryons in N to 
A transitions in Eq. (|14p . They are built from quark 
field operators with appropriate quantum numbers. For 
the nucleon (a spin-1/2 object) it is not unique. We 
consider a linear combination of the two standard local 
interpolating fields for the baryon octet: 

rjPiuud) = -2e'''"'[{u''^Cj5d'')u'' + P{u°''^Cd'')j5u''], 
r]"{ddu) = -2e'''"'[{d''^Cj5u^)d'' + l3{d''^Cu'')-f5d% 

r]^{uds) = -2y^e'^'«=[2(M'^^C75d'')s'= + K^C75s'')d'= 

-(d'^^C75s'')M" /3(2(w°'^Cd^)75s'= 
-H(u'^^Cs'')75d" - (d''^Cs'')75u'=)], 
ri^^iuus) = -2e'^''=[(M"'^C75s'')u'= + /3K^Cs^)75U=], 
ri^^uds) = -V2e'^''^[(u"'^C75s'')d=-h(ci"^C75s'')M'= 

+/3((u'^^Cs'')75d'= + (d"'^Cs'')75U^)], 
T]^'{dds) = -2e"''=[((i"^C'75S^)d'=-f/5(d°^ (75^)75^1, 
r]^\ssu) = -2e"''"[(s'^^C75u'')s'=-h/3(s''^CM^)75s'=], 
r]^~{ssd) = -2e"''"[(s"^C75d^)s'=-f /3(s°^Cd'')75s'=]. 

(18) 



Va i^ud) = ^e°'"=[2(M°^C7ad^)M'=-t- (u''^C7„u^)d'=], 
vf{ddu) = ^e°''"[2(d''^C7au'')d'= + (d''^C7„d'')M"], 
rj^'^iuus) = ^e°'"=[2(u°^C7„s^)u'= + (u°^C7„M^)s'=], 
T]^'\uds) = ^e°'"=[(w°'^C7ad'')s^-H (d°^C7as'')it'= 

+ {s-^Cj^u'')d% 
r]^' (dds) = ^e°'"=[2(d"^C7aS^)d'=-f (d'^^C7„d'')s^], 

Tjf\ssu) = ^e°^^[2(s''^C7„u^)s^ + (s"^C7aS^)u^], 
Vf'issd) = ■j^e^'>^[2{s-^Cja.d')s^ + is^^C-fc,s'>)d% 

(19) 



I 

B. Phenomenological representation 



We begin with the structure of the correlation hmction in the presence of the electromagnetic vertex to first order: 
n^ip) = 1 1 d'xe^P^ < 0\v^{x)H J d^yA,{y)f{y)]f]''mO > ■ (20) 
After inserting two complete sets of physical intermediate states, it becomes 

y.A^{y) < 0|?7^(a;)|Afc's' >< Ak' s'\j''{y)\Nks >< Nks\fi{0)\0 > . 
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Using the translation invariance on r]^{x) and j^{y), and the overlap strengths (Aa and Aat) of the interpolating fields 
with the states defined by, 

<0|7y„(0)|AA:',s' >= \^u„{k\ s'), ,^2) 
< 7Vfcs|?7(0)|0 >= A^u(fc,s), ^ ' 

we can write out explicitly the contribution of the lowest lying state and designate by ESC the excited state contri- 
butions: 

Af,{y)e'^P-'''y-'e"^y J2 M^' , s')up{k' , sOr"^ E u{k, s)u{k, s) + ESC. ^^^^ 

s' s 

The spinor sum for N in the above expression is 

u{k, s)u{k, s) = k + rriN, (24) 

and the spinor sum for A is 

^ w„(fc', s')up{k', s') ^~{k' + mA){9ap ~ llalp ^ + ^^^^^^i' (25) 

s' ^ 

We work in the fixed-point gauge in which A^{y) — —\F^^y'^. After changing variables from k to q = k' ^ k, which 
leads to d'^k ~ —d'^q, we have 



n„(p) = -X*^X^F^, J d^xd^y^^W^ - ml][{q - k'f - m%]-^ 

xe^(p-'c>(-l^e^«y)(fc' + mA){ggp - ^lalp - + %;^} x T'^'^ x (fc' - q + m^). 



(26) 



Integrating by parts and using properties of S functions, the expression collapses to 
n„(p) = l^F,4p^-ml]-^[p^~m%]-\p + mA) 



2 

X [Gi{gp,yjp. - j^gpphs + G2{gpuPp - Pvgpp)lb\{p + m^v), 



x{.9„.-i7"7p-li--^^-™-} (27) 



where we have chosen to work with the transition amplitudes Gi and G2 as defined in Eq. ^ and Eq. ([3]), rather 
than Gmi and Ge2- They are readily converted back and forth by Eq. ([7]) and Eq. (|5]). 

A straightforward evaluation of this expression leads to numerous tensor structures, but not all of them are inde- 
pendent of each other. The dependencies can be removed by ordering the gamma matrices in a specific order. After 
a lengthy calculation we identified 12 independent tensor structures. They can be organized (aside from an overall 
factor of i that also appears on the QCD side) as: 

n„(p) = _ ™^]-i[p2 _ 

[WEidGimA+fGimw) 

+WE2(|p2GimXi-4/3GiTOA-fGiTOjv) 

+WE3(fGimA+|Gimw) 

-HWE4(-|p2G'iTOX^-F|p2G2-h|GimA+|GiTOjv-|G2TOAr7ijv) 
-hWE5(-|GiTOX^mA,-|G2TOX^mAr-h|G2) 

+WE6(2p'G2-2G2TOATOjv) (28) 

+W0i(-|p2Gi-|GimAmAr) 

-HW02(|GiTOX^m7v+|Gi) 

+W03(-|p2Gi-|GimAmjv) 

+W04(|p'Gim^2^|p2G2mX'+|GimX'mAr-|Gi-|G2TOA+fG2mAr) 

+W05(-|GiTOX^TOAr-2Gi-|G2mA+|G2mjv) 

-hW06(lGirnX^TOAr-fG2TOA+|G2mAr)]. 
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The tensor structures associated with WEi and WOi are listed as follows: 

WEi = F^^ulbPltilwla WOi = F^i/ 757^1 71^70 

WE2 = Ff_,^-f5-ff^J^Pa WO2 = Ff^^j5pjf_,-ft,Pa 

WE3 = F^.^-y^pi'y^gaf, - J^gau) WO3 = -F^i.75(7i'SaM - 7m5qi/) (-29) 
VF£'4 = F,i^75(7^7ap^ - 7i,7qP^) = F^^j^ij^p^ - 7^p^)pa 

WE'S = Fi_,^j5p{jf_,p^Pa - ^„P^Pa) WO5 = F^,,j5p{j^JaP^ - ■^ulaPfi) 
WEq = F^_,^-fz{ga^_,Pu - gauP^i) WOe = Ff^^y"fzp[ga^,Py - gauPfi) 

The naming convention of M^E'i and VFOi is motivated by even-odd considerations. It turns out that the QCD sum 
rules associated with WE.^ will have even-dimension vacuum condensates, while the ones associated with WOi will 
have odd-dimension ones. Interestingly, the WEi structures themselves contain odd number of gamma matrices, while 
WOi even number. 

The next step is to perform the standard Borel transform defined by 

hm l(V)"+^(;^)'7(/), (30) 

~p .n — *oo il" (^P 

which turns a function of p^ to a function of A'P {M is the Borel mass). The transform is needed to suppress 
excited-state contributions relative to the ground-state. Unlike the N to N transition in the calculation of magnetic 
moments [lo| . the double-pole structure in the N to A transition in Eq. (j28p has unequal masses ttt-a and itln. The 
existence of an equal-mass double pole in the ground state was crucial in isolating the magnetic moments from the 
excited states because it gives a different functional dependence on the Borel mass (e-™' /*^' /A/2) than the single 
poles (e~™ ). Here for N to A transitions, we shall approximate the unequal-mass double pole by an equal-mass 
double pole, 

11 1 

o)', (31) 



p2 — TOAT p2 — mA p^ ~ ii- 

where = (m^ -I- m^)/2 can be considered as an averaged square mass. For the physical masses involved in the 
N to A transitions, this approximation is good to within 3% in the Borel region of interest we are going to explore. 
Because of this approximation, the determination of the N to A transition amplitudes will not be as accurate as 
the magnetic moments in the QCD sum rule approach. With the introduction of the pure double pole, the Borel 
transform of the ground state in Eq. (pS)) is given by 

[WEifGi(mA+mN) 

+ WE2lGi^{~NP + m^- 2m\ - rriNm^) 
+ WE-ilGi{mi^+mN) 

- WEi{lGi^{-M^ + m2 - Am\ - niNm^) + f G2(-M2 + - mA^jv)) 

- WE^i^G^^m^ + |G2^(-m^ + m^)) 

+ WEe2G2{-M^ + -niAmN) (32) 

- WOi^Gii-AP + fh'^ + mAniN) 
+ W02lGi^{mN + mA) 

- WOs^Gii-lVP + m2 + mAmw) 

+ W04{IGi^{-AP + fh^+ niNniA - ml) + IG2^{~AP + rh^ ~ 2m\ + mjvmA)) 

- WOz{lGi-;^(mN + tua) - |G2(toa - ^Tv)) 
3^ 



WOeilGij^mN - |G2(mA - mjv))], 



where M is the Borel mass parameter, not to be confused with the particle masses m^y, mA, or m. 
The excited states must be treated with care. The pole structure can be written in the generic form 

Ci I C2 I C3 I d ^_ /qq\ 

(p2_m^)(p2_„2^) Z.^ (p2-m^.)(p2_„i2^) (p2-m^)(p2_„^j ^^a. )(P^-'"i. ) ' (33) 
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s ^ / f 

/ / / / 

/ / / / 

Z + Z + Z + Z 

FIG. 1: The four contributions to the spectral function in the presence of an external field: transition between ground state N 
to ground state A, transitions between ground state N and excited states A*, transition between excited state A'^* and ground 
state A and pure excited states A''* A*. 



where Ci are constants. The first term is the ground state pole which contains the desired amplitudes G\ and Gi . The 
second and third terms represent the non-diagonal transitions between the ground state and the excited states caused 
by the external field. The forth term is pure excited-state to excited-state transitions. These different contributions 
can be represented by the diagrams in Fig. [1] Upon Borel transform, it takes the form 



C4 „-rn'^/M^ 



+e 



C2{l-e- 



A* 



C3(l- 



(34) 



The important point is that the non-diagonal transitions 
{N* to A and A^ to A*) give rise to a contribution that is 
not exponentially suppressed relative to the ground state. 
This is a general feature of the external-field technique. 
The strength of such transitions at each structure is a pri- 
ori unknown and is a potential source of contamination in 
the determination of the amplitudes. The standard treat- 
ment of the off-digonal transitions is to approximate it 
(the quantity in the square brackets) by a parameter A, 
which is to be extracted from the sum rule along with the 
ground state property of interest. Inclusion of such con- 
tributions is necessary for the correct extraction of the 
amplitudes. The pure excited state transitions (A^* to 
A*) are exponentially suppressed relative to the ground 
state and can be modeled in the usual way by introducing 
a continuum model and threshold parameter w. 

In summary, we identified 12 different tensor structures 
from which 12 sum rules can be constructed, once the 
correlation function is evaluated at the quark level. 



(with uud quark content) is: 

< 0\T{i^^^ {x)fjP {0)}\0 >= __|eahcga'&'='{ 
^r'TsC^f "C7a5-' + SfTri^^CS^^'Cj^Sf) 

+Sfj5Tr{CS™'C^c.Sf) - SfCS^'C^cS^/]}. 

(35) 

It is not necessary to list all the master formulas for the 
other transitions separately. Only three are distinct; oth- 
ers can be obtained by making appropriate substitutions 
in the above expression as specified below: 

• for ri7 — !■ A'' (with ddu quark content), interchange 
d quark and u quark in A"*" , 

• for S+7 — > E*+ (with uus quark content), replace 
d quark by s quark in p7 ^ A"'' , 



C. Calculation of the QCD side 

We start by writing down the master formulas used 
in the calculation of the QCD side, which are functions 
of the quark propagators, obtained by contracting out 
pairs of time-ordered quark-field operators in the two- 
point correlation function in Eq. (jlSp using the inter- 
polating fields in Eq. HH]) and Eq. We note that 
the same master formulas enter lattice QCD calculations 
where the fully-interacting quark propagators are gen- 
erated numerically, instead of the analytical ones used 
here. The master formula for the transition jp A+ 



• for S 7 — > S* (with dds quark content), replace 
u quark by d quark in S+7 —>■ 

• for S°7 — > (with ssu quark content), inter- 
change u quark and s quark in E"'"7 — > I]*^", 



• for ^ 7 — > ^* (with ssd quark content), replace 
u quark by d quark in S°7 

Here the interchange of two quarks is achieved by simply 
switching the corresponding propagators. 



The master formula for E''7 — > S*^ (with uds quark 
content) is: 

< 0|T{r;f" (:e)77^"(0)}|0 >= ^e^b-e'^'b'-'l 
-S^^'j^CS^'Cj^Sf ~ SfTr{^^CS™'C-f^SX') 

(36) 

The master formula for A7 — s- (with uds quark con- 
tent) is: 

< 0\T{7]^'\x)fj^}\0 >= |e'^'''=e°''''^'{ 

-2SfTr{j,CS^''''Cj^Sf) + Sfj^CS^'C-fo^Sf 
-ST-f^CS^^'C-f^Sf - Sr'Tr{-f^CS™'C^^S'i,-') 

+2S:i^CS™'CjaSfj5 - 2SfTr{CS™'Cj^Sf)j5 
~SfTr{CSr''Cjo.S^/h, - SfCSr'^'Cj^S^/j5 

(37) 

In the above equations, 

Sf{x, 0; F) = (0 I T{ q'^ix) <f{0) } \ 0)f, q = u,d, s, 

(38) 

is the fully interacting quark propagator in the presence 
of the electromagnetic field. It is obtained by the oper- 
ator product expansion (OPE) and is made of dozens of 
analytic terms. It has been discussed to various degree in 
the literature 0, [l^, [13, [2^ so we will not list the terms 
here. 

In addition to the standard vacuum condensates, the 
vacuum susceptibilities induced by the external field are 
defined by 

{q9cGf,^q)F = eqK{qq)Ff,iy, (39) 
{q9c^t,upxGP^-i^q)F = ieq^{qq)Ff,^. 

Note that x has the dimension of GeV^^, while k and ^ 
are dimensionless. 

With the above elements in hand, it is straightforward 
to evaluate the correlation function by substituting the 
quark propagator into the various master formulas. We 
keep terms to first order in the external field and in the 
strange quark mass (u-quark and d-quark masses are ig- 
nored). Vacuum condensates up to dimension eight are 
retained in the expansion. The algebra is extremely te- 
dious. Each term in the master formula is a product 
of three copies of the quark propagator, and there are 
hundreds of such terms over various color permutations. 
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FIG. 2: Diagrams considered for the calculation of N to A 
transitions that do not contain strange quark mass terms. 

The calculation can be organized by the diagrams (sim- 
ilar to Feynmann diagrams) in Fig. [5] and Fig. [31 Note 
that each diagram is only generic and all possible color 
permutations are implied. We used the computer alge- 
bra package called REDUCE to carry out some of the 
calculations. The QCD side has the same tensor struc- 
ture as the phenomenological side and the results can be 
collected according to the same 12 independent tensor 
structures in Eq. ((29)) . 

D. The QCD sum rules 

Once we have both the QCD side (commonly referred 
as LHS) and the phenomenological side (RHS), we can 
isolate the sum rules by matching the two sides. Since 
there are 12 independent tensor structures, 12 sum rules 
can be constructed for each transition. For the 8 transi- 
tions, the total number of sum rules is 96. We divide the 
sum rules into four groups according to their dependence 
on the amplitudes and performance. The first group con- 
sists of the sum rules at structures WEi, WO2 and WE3 
which only contain the amplitude Gi , and at WEg which 
only contains G2- These sum rules are expected to per- 
form better because they involve fewer parameters. The 
second group consists of the sum rules at structures WE5, 
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FIG. 3: Diagrams considered for the calculation of N to A 
transitions that are proportional to the strange quark mass 
(denoted by x). 



WO5 and WOe which have both Gi and G2. It is diffi- 
cult to extract both amplitudes at the same time. One 
way out is to fix one of them and treat the other as 
free parameter. This reduces the predictive ability of the 
QCD sum rules. Our numerical analysis suggests that it 
can serve as a consistence check at best. Another way is 
to take appropriate combinations of sum rules to elim- 
inate one of the amplitudes. We tried but found the 
sum rules obtained this way perform worse than the un- 
manipulated ones. We advocate against the practice of 
manipulating the original sum rules like algebraic equa- 
tions (such as taking derivatives or linear combinations) 
to obtain new sum rules. The third group consists of the 
sum rules at structures WE2, WOi and WO3 which have 
one of the amplitudes and additional dependence on the 
averaged square mass rn^, aside from the common factor 
g-TO /M _ -y^g found these sum rules perform consistently 
worse than those in groups 1 and 2, possibly because of 
the approximation introduced by fh?. The fourth group 
consists of the sum rules at WE4 and WO4 which contain 
both Gi and G2, and additional dependence. They 
perform consistently the worst. In the following, we are 
going to give detailed analysis of the sum rules in group 
1, but first let us present the sum rules in this group. 

At the structure WEi , the sum rules can be expressed 



in the following form: 

ciL-4/27£:2(«;)M4 + C2xamsL-^"/^'^Ei{w)M'^ 
+C3amsL-^/^^Eo{w) + C4xa^L^/^^ Eo{w) 
+C5bL-^/^'Eoiw) + cemlam.L-^^/^^^ (40) 

+Cra^L^0/27_l^ + ,^^^2^2 j^- 10/27 

= 'Xn'Xab-"^"/^" [^GiiruA + tun) + A]. 

Here A at and A a are the rescaled current coupling 
A^ = (27r)^A^. The quark condensate, gluon condensate, 
and the mixed condensate are represented by 

a = — (2-77)^ (uu), b = {glG^), tuq = —{ugccr-Gu)/{uu). 

(41) 

The quark charge factors are given in units of electric 
charge 



e„ = 2/3, ed = -1/3, = -1/3. 



(42) 



Note that we choose to keep the quark charge factors 
explicit in the sum rules. The advantage is that it can 
facilitate the study of individual quark contributions to 
the amplitudes. The parameters / and account for 
the flavor-symmetry breaking of the strange quark in the 
condensates and susceptibilities: 



/ 



(ss) ^ {s(j,,a ■ Gs) j^^Xs_ 
{uu) {uQca ■ Gu) ' X 



Kg 

K 
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(43) 



The anomalous dimension corrections of the interpolat- 
ing fields and the various operators are taken into account 
in the leading logarithmic approximation via the factor 



ln(MVA2 



QCD 



(44) 



where /x = 500 MeV is the renormalization scale and 
^QCD is the QCD scale parameter. As usual, the pure 
excited state contributions are modeled using terms on 
the OPE side surviving — > 00 under the assumption 
of duality, and are represented by the factors 



En{w) 



n! 



(45) 



where w is an effective continuum threshold and it is in 
principle different for different sum rules and we will treat 
it as a free parameter in the analysis. 

The Ci coefficients differ from transition to transition. 
It is not necessary to list the coefRcients separately for all 
8 transition channels. We only need to carry out three 
separate calculations for S+7 S*"*", 'EP'j S**^ and 
A7 S**^. Other channels can be obtained from them 
by making appropriate substitutions as specified below: 

• for p7 — > A+, replace s quark by d quark in S+7 

S*+, 

• for n7 ^ A", interchange d quark and u quark in 

proton p7 A+ , 



10 



• for E-7 ^ S* 



replace u quark by d quark in and in the transition channel A7 — > E*", 



• for add corresponding coefficients of 

S+7 ^ E*+ and E-7 ^ E*", then divide by 2, 



• for ^ 7 ^ S* , replace u quark by d quark in 



Here the interchange between u and d quarks is achieved 
by simply switching their charge factors e„ and e^. The 
conversions from s quark to u or d quarks involve set- 
ting m-s = 0, f — (f> — 1, in addition to switching the 
charge factors. In practice, we computed all 8 channels 
separately, and used the substitutions as a check of our 
calculations. 

The coefficients appearing in Eq. ( I40p are given by, in 
the transition channel E+7 E*+, 



ci = i(-l + /3)(e, -e„), 

C2 = -i(l + 2/3)(e„-eJ,0), 

C3 = i^(e„(-l + P + 3fs + Wis + K - /3k - (2 + m 

+e,(-6 + /,(4 -KC^ + 2^0 + /3(-6 + /s(2 + 0(« + 0)))) 
C4 = i(-eu(l + + fs) + (2 + P)esfs^), 
C5 = -^{-1 + P){es -e„), 
C6 = ^(l + 2/3)(2e,-e„(H-/3)), 

C7 = M(-e«(2 - + e + /3(2 + + + /s(-4 + + 0) 

+e«(-6 + /.(4 + 2/3 - + /3k0 + (2 + /3)0O)), 
C8 = n4(eu(8 + 2/3 + 8/3 + 5/3/,) - (16 + 7/3)e,/,</.); 

(46) 

in the transition channel S''7 S*°, 



cl = i(-(2 + I3)ed + e, + (2 + /3)e„), 
C2 = i((3-2/3)ed + 2(-l + /3)e,0, 
C3 = M(ed(-2 - 12/,s - At + /3(1 + 3/, - K - + 5^ 
+e„(-2 + 6/, + 5k - 4e + /3(-l - 3/, + k + 0) 

+e.(-6 + /,(-2 + 2K0-</)O)), 
C4- i^(-eu(2-3/3 + 2(l + /3)/,) 

+ed(3 + 2fs + /3(-3 + 2/,)) ~ Csfscb), 
C5 = Tfe((5 + 2/3)e<i + 3e, - (3 + 2/3)e„), 
C6 = |^("2e,, - (-9 + 4/3)erf(l + /,) + (-5 + 4/3)e„(l + /,)), 
C7 = ^(e./.(-2 + 2k(/. - (/-O + ed(2/,(-7 - 2k + 

+3(2 + K + + /5(-3 + 2/,)(2 + K + e)) - e„(2(2 + k + 

+/3(-3 + 2/,,)(2 + K + + fsi-S -7k + 20)), 
C8 = 2S8 (-f'rf(23 + 17/s + /3(-20 + 13/,)) 

+e„(15(l + /,) + /3(-20 + 13/,))). 

(48) 

At structure WEg, the sum rule is: 



+C3xa^L^/^'Eo{w) + C4a2i20/27_i^ + ^^^^2^2^-10/27 _i 

. . (^^^ 

The coefficients appearing in Eq.f l49p are given by, in the 
transition channel E+7 — )■ E*+, 



/3(-e„ + e,/,0), 

i/3(e„(-4 - 2k + + e,/,(4 + 2k(/. - (/-O), 



-|/3(e„ - e,/,(/)). 



Cl 

C2 = 
C3 = 

C4 = j^(_3e„(-l + /,)e + 2/3(e„(-4 - 2k + 

+e,/,(4 + 2K0-<^O)), 
C5 = g/3(e„ - e,/,0); 

in the transition channel S°7 ^ S**^, 



(50) 



Cl = -|(-l + /3)(e, -e„), 
C2 = |(1 + 3/3) (e„ - e,/,^), 

C3 = ^(e./. - e,/,0(K + e) + e„(2 - 3/, + k + C) 

+/3(e„(8 - 3/, + K + C) - e,(6 + /,(-l + (/)(k + e))))), 
C4 = i(l + /3)/,(e„ - es/s(/)), 

C5 = li2(~l + -e„), 

C6 = 3^(e,(2-/3(-4 + /,)-5/,) 
+e„(l + 2/,+/3(-7 + 4/,))), 

C7 = M(e«(~3(l + /3)/2 + 3(1 + /3)k 
+/,(2-2K + e + /3(8-2K + e))) 
+e,(-3 + /,(6 - 3k0 + /,(-2 + 2Kcf> - 0^)) 
+/3(-3 - 3/,K0 + /2(-2 + 2k0 - cbO))), 

cs = -134(1 + /3)(e„(-l + 8/,) - e,/,(2 + 5/,)</)); 

(47) 



Cl = /3(e„ - es/,(/)), 

C2 = i/3(e„(4 + 2k - + e,/,(-4 - 2k(/. + 0^)), 
C3 = |/3/s(e„ - e,/,(/)), 

C4 = j^/,(3e,(-l + /,)</>e + 2/3(e„(4 + 2k - 

+e,/,(-4~2K,/) + <^0)), 
C5 = g/3/s(-e„ + es/s0); 

and in the transition channel A7 —> E*'', 

Cl = i/3(ed - e„), 

C2 = i/3(ed-e„)(4 + 2K-e), 

Cs = i/3(ed - e„), (52) 
C4 = jgied - e„)(/3(8 + 4k - 2^ - 3(-l + /,)0, 
C5 = i^/3(-e<i + e„). 
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At structure WO2, the sum rule is: 



and in the transition channel A7 — > E 



(53) 



The coefficients appearing in Eq. ( I53p are given by, in 
the transition channel E+7 



C2 = ^(-e„ + ejs4>), 

C3 = i^(e„(-2 + 6/, + 2k + 6/3(-2 + 2/, + n) - ^ 

+e,(-6 + /,(2 - 2(1 + 3/3)/t</) + </.C))), 
C4 = j^(4e, - e„(l + 2/3(-l + /,) + 3/,)), 
C5 = + + /s) + 3/s) - 4e 

C7 = 3^(e,(6 + /,(-4 + (/.(k + 6/3k - 20)) 

+e„(2 - 2k + e + /s(-4 + AC - 6/3« + 0)), 
C8 = ii4(-2e« + e,(l - 2/3(-l + /,) + /,)); 



(54) 



ci = i(l+/3)(e, -e„), 

C2 = - i^(-l + ^)(e« - es/s(/)), 

C3 = ^(e„(2 - 6/, + K + - e./.(-4 + ,^(k + 0) 

+/3(e,(-14 + 6/3-7K-0 

+e,(24 + /,(-16 + 7Ac0 + 0e)))), 



^4 = ^(-2(es + QPes 



6l3e^) 



cl = |(-3ed + Bs + 2e„), 
C2 = ^(ed - e„), 

C3 = ^{6es + e„(2 + 6/, + 7k- 6/3(-2 + 2/, + k) - 2^ 
+erf(-8 - 6/, - 7k + 6/3(-2 + 2/, + k) + 2^), 



Ucdi^O + P{23 - 22fs) + 16fs 
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+e„(-2(7 + 8/,) + /5(-23 + 26/,))), 
C5 = i(e<i(-3 + 2/3(-l + /.)-2/,) 
+2e„(l + /3 + /,-/3/,) + e,/,(/.), 

C6= 2i8((-l + 2/3)erf-(l + 2/3)e«), 
C7 = ^(-6es - e„(K - 12/3k 

+/,(8+ (4 + 6/3)^-20-2(2 + 0) 

+ed(/.(14 + K + 6/3k - 20 - 2(2 + (-2 + 6/3)k + 0)), 
C8 = T]k(l - 8e,, + e„(6 + /3(7 - 10/,) + 8/,) 

+ed{~U-7p-8fs+6(3fs)). 

(56) 

This concludes the presentation of the QCD sum rules 
that will be used in the numerical analysis. 

IV. SUM RULE ANALYSIS 

A. General procedure 

The sum rules for N to A transition have the generic 
form of OPE - ESC = Pole + Transition, or 



Utran{QCD,/3,W,AP) = 



(57) 



+ ((-5 + 7/3)e, + (9-7/3)e„)/,), 
C5 = ^(j3{-'^ + fs){2eu-7esfscf>) 

+/s(-4e„ + e,(l + 3/,)0)), 
C6 = -218 + 3/3)(e„ - esfs<l>), 

C7 = ^(e,(12 + 2/, (-8 - 5/3 + /, + (2 + /3(-5 + /,) - /,) 
K(t>) - (-1 + /3)/s(l + /.)</'0 + 2e„(-3K + fs{-2 + 3/, 
+2k - + /3(2 + k + 3/,(l + k) + 0)), 

C8= 5ig(-2(-l + 2/3)(e, +e„) 
+ (e. + 5/3e, + (-5 + 3/3)e„)/,); 

(55) 



where QCD represents all the QCD input parameters. 
/(Gi, G2) is a function of transition amplitudes Gi and 
G2. The basic mathematical task is: given the function 
Iltran with known QCD input parameters, find the phe- 
nomenological parameters (Gi, G2, transition strength 
A, coupling strength A, and continuum threshold w) by 
matching the two sides over some region in the Borel mass 
M. A minimization is best suited for this purpose. It 
turns out that there are too many fit parameters for this 
procedure to be successful in general. To alleviate the 
situation, we employ the corresponding mass sum rules 
which have a generic form of OPE - ESC = Pole, or 

Il^assMQCD,(3,w,,M^) = ~X%e-"'-/^'"; (58) 
U^ass,AiQCD,P,W2,M^) = Aie-"A/M^. (59) 

They share some of the common parameters and factors 
with the transition sum rules. Note that the continuum 
thresholds may not be the same in different sum rules. 
By taking the following combination of the transition and 
mass sum rules. 



,^{QCD,l3,w,M^) 



y/U„,ass.NiQCD,0,wuM^)n„,ass &(QCD,p,W2.M^) (60) 
_ /(Gi G2) ,4 
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the couplings A and the exponential factors are canceled 
out. This is another advantage for introducing the ifi^ 
prescription. The form in Eq. ([60|) is what we are going 
to implement. By plotting the two sides as a function of 
1/M^, the slope will be related to the transition ampli- 
tudes Gi or G2 and the intercept the transition strength 
A. The linearity (or deviation from it) of the left-hand 
side gives an indication of OPE convergence and the role 
of excited states. The two sides are expected to match 
for a good sum rule over a certain window. This way of 
matching the sum rules has two advantages. First, the 
slope, which is proportional to the transition amplitudes, 
is usually better determined than the intercept. Second, 
by allowing the possibility of different continuum thresh- 
olds, we ensure that both sum rules stay in their valid 
regimes. 

For the octet baryons, we use the chiral-even mass sum 
rules in Ref. [23| which are reproduced here in the same 
notation used in this work. 

+P4a^L^/^ + p^a'KL^'^ + p^mla'L-^/^^ ^ 

(61) 

The coefficients are, for N: 

Pi = ^(5 + 2/3 -1-5/32), p2-_i_(5 + + 5/32), 
P3 = 0,p4 = ^(7-2/3-5/3), P5 = 0, 
-^(13 -2/3 -11/32); 

(62) 

for A: 

pi = ^(5 + 2/3 + 5/32), 
P2 = 556 (5 + 2/3 + 5/32), 

Pz = ^((20 - 15/,) - (16 + 6/,)/3 - (4 + 15/,)/32), 
Pi = - 5 - 6<) + (4 + 4/,)/3 + (4/, + 1 + 6t)/32), 

P5 = ^((10/s + 11) + (2 - 8/,)/3 - (2/, + 13)/?2), 
P6 - 2S8 ((-16/s - 23) + (8/, - 2)/3 + (8/, + 25)/32); 

(63) 

for S: 

pi = 5^(5 + 2/3 + 5/32), 

P2 = 2^5(5 + 2/3 + 5/32), 

P3 = 3l((12 - 5/.) - 2/,/3 - (12 + 5/,)/32), 

PA = -^((4/s + 21 + 18<) + 4/,/3 + (4/, - 21 - 18i)/32), 

P5 = 2l((6/s + l)-2/3-(6/,-l)/32), 

P6 = -m((12./s + 1) - 2/3 - (12/, - l)/32); 

(64) 

for 

PI = 5^(5 + 2/3 + 5/32), 

P2 = 2i6(5 + 2/? + 5/32), 

P3 = i|((2 - /.) - 2/,/3 - (2 + /,)/32), 

Pi = -^((15 - /. + m) - 10/,/3 - (15 + /, + 18t)/32), 

P5 = j-JsiUs + 6) - 2/,/3 + (A - 6)/32), 

P6 - -^Jsiifs + 12) - 2/,/3 + (/, - 12)/32). 

(65) 



The function t is defined as t = In ^ —Jem with jem ~ 
0.577 the Euler-Mascheroni constant. 

For the decuplet baryons, we use the chiral-odd sum 
rules [25! at the structure g^^, which are reproduced here, 
for A: 

I a El ii6/27 Mi - I mla Eq L2/27 M^~^ab 2.16/27 

(66) 

for E*: 
|(/, + 2)ai?iLi6/2VAf4 

-§(/, + 2) m2a i?o i2/27 M2 - ^ (/, + 2) a6 L^e/^T 

+ i TO, E2 L-8/27 _^ 2 ^^^2 ^16/27 

= A2. M,. e-^'^^-/*^'; 
for E*: 

I (2/, + 1) a i?i m4 

-I (2/, + 1) TO2a Eo ^2/27 M2 - ^ (2/, + 1) a 6 LI6/27 

+TO, E2 L-«/2^ Af6 + I TO, /, K,a2 ^16/27 

= A2. M,. e-^'^B*/A^^ 

B. QCD input parameters 

The standard vacuum condensates are taken as a = 
0.52 GeV^, b = 1.2 GeV^, to§ = 0.72 GeV2. We take into 
account the possible factorization violation in the four- 
quark condensate (in terms of k^jO?) and use k„ — 2.0. 
The strange quark parameters are placed at to, = 0.15 
GeV, / = 0.83, (f) = 0.60 [H, [13, 123. The QCD scale 
parameter is restricted to Aqcd = 0.15 GeV. This set of 
parameters is used in all studies of hadron properties in 
the QCD sum rule approach. In the presence of exter- 
nal fields, additional condensates called vacuum suscep- 
tibilities are introduced. These parameters are less well- 
known. They have been estimated in studies of nucleon 
magnetic moments [1, [13, do, 19] where the availability 
of precise experimental data has put strict constraints on 
these parameters. In particular, the magnetic suscepti- 
bility X is the subject of two recent studies [30, 31]. We 
use the values x = —6.0 GeV~'^ and n = 0.75, ^ — —1.5. 
Note that x is almost an order of magnitude larger than 
K and ^, and is the most important of the three. 

The above parameters are just central values. We will 
explore sensitivity to these parameters by assigning un- 
certainties to them. To this end, we use the Monte-Carlo 
procedure first introduced in Ref. [3?] which allows the 
most rigorous error analysis of QCD rum rules. The 
procedure explores the entire phase-space of the input 
QCD parameters simultaneously, and maps it into un- 
certainties in the phenomenological parameters. It goes 
briefly as follows. First, a sample of randomly-selected, 
Gaussianly-distributed condensates are generated with 
assigned uncertainties. Here we give 10% for the un- 
certainties of input parameters, and this number can be 
adjusted to test the sensitivity to the QCD parameters. 
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Then the OPE is constructed in the Borel window with 
evenly distributed points Alj. Note that the uncertain- 
ties in the OPE are not uniform throughout the Borel 
window. They are larger at the lower end where uncer- 
tainties in the higher-dimensional condensates dominate. 
Thus, it is crucial that the appropriate weight is used 
in the calculation of ■ For the OPE obtained from the 
k'th set of QCD parameters, the per degree of freedom 
is 

xl _^ [nr-(M|)-nr" (M|)]^ 

Ndf fr{ {nB-np)al,,{M,) ' ^ ' 

where n°^^ refers to the LHS of Eq. §Q and W'"" 
its RHS. The integer Up is the number of phenomeno- 
logical search parameters. In this work, nB=5l points 
were used along the Borel window. The procedure is re- 
peated for many QCD parameter samples, resulting in 
distributions for phenomenological fit parameters, from 
which errors are derived. In practice, 200 samples are 
sufficient for getting stable uncertainties. We used about 
2000 samples to resolve more subtle correlations among 
the QCD parameters and the phenomenological fit pa- 
rameters. This means that each sum rule is fitted 2000 
times to arrive at the final results. 



V. RESULT AND DISCUSSION 

We did a wide survey of the 96 sum rules at the 12 
structures and the 8 transitions. Overall, based on the 
quality of the match, the broadness of the Borel window 
and its reach into the lower end, the size of the con- 
tinuum and non-diagonal contributions, and OPE con- 
vergence, we found that most of the sum rules do not 
perform well in terms of their ability to predict stable 
values for the amplitudes. They can be roughly divided 
into four groups using our naming convention of WE^ 
and WOi, as alluded to earlier. Group one consistently 
produces the most reliable results, except the structure 
at WE3. In the following, we focus on the sum rules at 
WEi and WO2, which only contain the amplitude Gi, 
and at WEg which only contains G2. 



A. The sum rules at WEi 

The sum rules are found in Eq. (PD)) . Ideally, we would 
like to extract 6 parameters: Gi, A, w, wi,W2 and [3. But 
a search treating all six parameters as free does not work 
because there is not enough information in the OPE. In 
fact, the freedom to vary (3 can be used as an advan- 
tage to yield the optimal match. One choice is (3 — —0.2 
which minimizes the perturbative term in the mass sum 
rule (the first term in Eq. jssj . Another parame- 

ter that can be used to our advantage is the continuum 
thresholds wi and 1112 for the corresponding mass sum 



rules. We fix them to the values that give the best solu- 
tion to the mass sum rules independently. The following 
values are used: for the nucleon, wi — 1.44 GeV; for A, 
wi = 1.60 GeV; for S, wi = 1.66 GeV; for S, wi = 1.82 
GeV; for the A, W2 = 1.65 GeV; for S*, u;2 = 1-80 GeV; 
for S*, u'2 = 2.0 GeV. In this way the transition and the 
mass sum rules can stay in their respective valid Borel 
regimes. This leaves us with three parameters: Gi, A, 
w. Unfortunately, a three-parameter search is either un- 
stable or returns values for w smaller than the particle 
masses, a clearly unphysical situation. Again we think 
this is a symptom of insufficient information in the OPE 
to resolve the parameters. So we are forced to fix the con- 
tinuum threshold w to the value that corresponds to the 
best match for the central values of the QCD parameters, 
and extract two parameters: Gi and A. 

TABLE II: Results for Gi in A ^ A transitions from the 
QCD sum rule at structure WEi in Eq. (|40)) . The columns 
correspond to, from left to right: transition channel, l3 value, 
Borel window, continuum threshold, transition strength A, 
amplitude Gi . The errors in the results are derived from 2000 
samples in the Monte-Carlo analysis with 10% uncertainty on 
all the QCD input parameters. 



Transition 


P 


Window 


w 


A 


Gi 








(GeV) 


(GeV) (GeV-2) 


(GeV-i) 


PI ' 


A+ 


-0.2 


0.95 to 1.3 


1.4 


-0.54(7) 


3.84 (37) 




AO 


-0.2 


0.95 to 1.3 


1.4 


-0.54(7) 


-3.84 (37) 


E+7 


-.E*+ 


-0.2 


1.2 to 1.4 


1.7 


0.07(1) 


2.74 (25) 


E«7 


E*° 


-0.2 


1.1 to 1.6 


1.6 


-0.02(1) 


1.18 (10) 


E-7 


^ E*- 


-0.2 


1.2 to 1.8 


1.9 


0.01(1) 


-0.33 (5) 


H«7 


^ 


-0.2 


1.2 to 1.8 


1.8 


0.04(1) 


-1.02(10) 


H~7 




-0.2 


1.2 to 2.0 


2.1 


0(1) 


0.10 (2) 


A7 - 


-^E*« 


-0.2 


1.0 to 1.4 


1.65 


-0.01(1) 


2.92(25) 



The results of such an analysis are given in Table [TTl 
The Borel window is determined by the following two cri- 
teria: OPE convergence which gives the lower bound, and 
ground-state dominance which gives the upper bound. It 
is done iteratively: using the optimal value of /3, we ad- 
just the Borel window until the best solution is found and 
the two criteria are roughly satisfied. We also checked 
that the results are not sensitive to small changes in (3 
and the Borel window. Our result on the A7 S*" tran- 
sition is a prediction: we know of no other calculations 
of this transition channel. The results on the parame- 
ter A indicate that the non-diagonal transitions are not 
significant in this sum rule, perhaps due to cancellations 
in the excited states, but they are not negligible in the 
proton and neutron channels. Ignoring them will alter 
the slope of the RHS and lead to different results for Gi . 
In fact, we found that one of the main symptoms of a 
sum rule performing poorly is the relatively large con- 
tribution of the non-diagonal transitions represented by 
A. The effects of such contributions are not known until 
a numerical analysis is carried out. We stress that the 
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errors are derived from Monte-Carlo distributions which 
give the most realistic estimation of the uncertainties. 
An example of such distributions is given in Fig. H) We 
see that they are roughly Gaussian distributions. The 
central value is taken as the average, and the error is 
one standard deviation of the distribution. We found 
about 10% accuracy in our Monte-Carlo analysis, result- 
ing from 10% uniform uncertainty in all the QCD input 
parameters. Of course, the uncertainties in the QCD pa- 
rameters can be non-uniform. For example, we tried the 
uncertainty assignments (which are quite conservative) 
in Ref. [32l |. and found about 30% uncertainties in our 
output. 




180 - 
160 - 
140 - 
120 - 
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FIG. 4: Histogram for the Gi transition amplitude (top) and 
transition amplitude (bottom) obtained from Monte-Carlo fits 
of Eq.( go} at WEi for 2000 QCD parameter sets. They 
are based on 10% uncertainty given to all the QCD input 
parameters. 

To gain a better appreciation on how the QCD sum 
rules produce the results, we show some details of the 
analysis in Fig.O using the transition p7 — > A+ as an ex- 
ample. There are three graphs in this figure to give three 
different aspects of the analysis. The first graph shows 
how the two sides of Eq. ([SO)) for this sum rule match over 
the Borel window, which should be linear as a function of 
1/M^ according to the right-hand side of this equation. 
Indeed, we observe good linear behavior from the LHS 
(which is OPE-ESC). The slope is directly proportional 
to the transition amplitude Gi, and the intercept give the 
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FIG. 5: Analysis of the QCD sum rule in Eq. (|40|l (structure 
WEi) for the proton at /3 = -0.2 according to Eq. (|60|- Top 
figure, the pole plus transition terms (solid lines) are com- 
pared against the OPE minus the excited-state contributions 
(dashed lines) as a function of 1/M^ (the two should match 
for an ideal sum rule) . Also plotted are the individual contri- 
butions from u (long-dashed lines) and d (dot-dashed lines) 
quarks. Middle figure, the total in the OPE side and its var- 
ious terms are plotted as a function of M^. Bottom figure, 
the 3 terms in the phenomenological side: pole (solid), tran- 
sition (long-dashed), and excited (dot-dashed) are plotted as 
a function of . 



non-diagonal transition contribution A. We also plot the 
individual contributions from u and d quarks. We see 
that in this transition, the u-quark contribution is the 
dominant one, which is expected because it is doubled- 
represented in the proton (uud). In the next two graphs, 
we give details on the individual terms in the sum rule in 
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Eq. (|40p , remembering that the sum rule is in the generic 
form of OPE-ESC=Pole+Transition. The second graph 
in Fig. [5] shows how the various terms contribute to the 
OPE as a function of AP. The M'^ terms, which contain 
the contributions from the condensates xa^ and 6, play 
an important role. It is followed by the term which 
is the perturbative contribution. The term is zero 
in this channel because it is proportional to the strange 
quark mass. The term gives a negative contribution 
and its contribution as a percentage of the entire OPE is 
about 10% at the lower end of the Borel window. In the 
third graph we show the three terms that comprise the 
phenomenological side (Pole, Transition, and ESC) as a 
function of M^. The ground-state pole is dominant at 
the low end of the Borel window. The excited-state con- 
tribution starts small, then grows with M^, as expected 
from the continuum model. The transition contribution 
is small in this sum rule. It is consistently smaller than 
the excited-state contribution and has a weak dependence 
on the Borel mass. 

In our Monte-Carlo analysis, the entire QCD input 
phase space is mapped into the phenomenological out- 
put space, so we can also look into correlations between 
any two parameters by way of scatter plots of the two pa- 
rameters of interest. Fig. [6] shows such an example. We 
see that the transition amplitude Gi has a strong nega- 
tive correlation with the vacuum susceptibility x- Larger 
X (in absolute terms since x is negative) leads to smaller 
Gi. Precise determination of the QCD parameters, es- 
pecially those that have strong correlations to the out- 
put parameters, is crucial for keeping the uncertainties in 
the spectral parameters under control. We found similar 
strong correlations with x in other transition channels. 

5.5 I . • • • • • • • 1 



5 - 




2.5' ' • • • • • • • ' 
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FIG. 6: A scatter plot showing correlations between the tran- 
sition amplitude Gi and one of the QCD parameters x for 
the P7 —* at structure WEi. It is obtained from 2000 
Monte-Carlo samples with 10% uncertainty on all the QCD 
parameters. 



B. The sum rule at WO2 

The sum rule at this structure is found in Eq. (|53|) . 
We use the same procedure to analyze it. Since OPE 
expansion in this sum rule goes deeper than others (from 
to l/M"*), it is expected to be more reliable than 
the WOi sum rule. But our analysis shows that this 
advantage is offset by the smallness of the 1/M^ and 
terms. As a result, its performance is about the 
same as the sum rule at WE 1. Table [ml displays the 
results extracted from this sum rule. 



TABLE III; Similar to Table HH but for the QCD sum rule in 
Eq. ([531) (structure WO2). 



Transition 


/3 


Window 


w 


A 


Gi 






(GeV) 


(GeV) (CeV-3) (GeV-^) 


P7 A"*" 


-0.2 


0.9 to 1.2 


1.3 


-0.5(1) 


3.52(32) 


nj ^ A° 


-0.2 


0.9 to 1.2 


1.35 


0.5(1) 


-3.68(33) 


E+7 E*+ 


-0.2 


1.2 to 1.6 


1.4 


0(1) 


2.20(20) 


E«7 S*« 


-0.2 


1.0 to 1.8 


1.6 


0.02(1) 


0.92(10) 


E-7 ^ E*- 


-0.2 


1.0 to 1.6 


1.6 


-0.01(1) 


-0.32(5) 




-0.2 


1.0 to 1.8 


1.65 


0.01(1) 


-0.91(10) 




-0.2 


1.2 to 1.8 


1.8 


0(1) 


0.21(3) 


A7 ^ E*" 


-0.2 


1.0 to 1.6 


1.45 


0.05(1) 


0.89(10) 



C. The sum rule at WEe 

The sum rule is found in Eq. This is the only sum 
rule that depends on G2 alone. The results of our analysis 
are given in Table [IVj The off-diagonal contributions (A) 
in the p — > A+ and n — > A° are smaller by half than those 
for Gi, but the uncertainties are larger. They are small 
in the other channels, just like those for Gi. rediction, 
but others are not good. The reason of this is because he 
transition contribution (A) are as larger as the G2 term, 
so the two terms can't distinguish and get contaminated. 



TABLE IV: Similar to Table [iVl but for the G2 amplitude 
from the QCD sum rule in Eq. ([49)l (structure WEe). 



Transition 


/3 


Window 


w 


A 


G2 








(GeV) 


(GeV) (GeV-^) 


(GeV-2) 


PI 


^ A+ 


-0.2 


1.0 to 1.3 


1.3 


-0.26(4) 


-2.34(48) 


717 


^ A" 


-0.2 


1.0 to 1.3 


1.3 


0.26(4) 


2.34(48) 


E+7 


^E*+ 


-0.2 


1.4 to 1.6 


1.8 


-0.07(1) 


-0.41(45) 


E«7 


^E*« 


-0.2 


1.5 to 1.8 


1.5 


-0.01(1) 


-0.31(16) 


E-7 


E*- 


-0.2 


1.2 to 1.6 


1.8 


0.01(1) 


-0.32(12) 


H«7 


^ w*o 


-0.2 


1.5 to 1.8 


1.8 


0.03(1) 


0.38(29) 


=-7 




-0.2 


1.4 to 1.8 


1.9 


0(1) 


0.18(7) 


A7 


-^E*° 


-0.2 


1.0 to 1.3 


1.45 


-0.09(1) 


-0.62(17) 
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TABLE V: Individual quark contributions to the transition 
amplitude Gi in units of GeV~^ extracted from the QCD 
sum rules in Eq. (140 |l (structure WEi). 









/-id 


Lri 




PI - 


A + 


2.00(25) 


1.28(13 ) 





3.84(38) 


nj 


^ AO 


-1.28(13) 


-2.56(25) 





-3.84(38) 


E+7 


->E*+ 


2.06(17) 





0.69(6) 


2.75(22) 


S«7 




1.00(9) 


-0.50(5) 


0.67(6) 


1.17(10) 


E-7 


-+ s*- 





-1.04(8) 


0.71(6) 


-0.33(6) 


=07 




-0.71(6) 





-0.25(2) 


-0.97(7) 


="7 







0.37(3) 


-0.27(2) 


0.11(2) 


A7 - 




1.91(16) 


1.17(10) 


-0.16(2) 


2.92(25) 



TABLE VL Individual quark contributions to the transition 
amplitude G2 in units of GeV~^ extracted from the QCD sum 
rules in Eq. (|49)l (structure WEe). 















VI - 


A + 


-1.57(32) 


-0.77(10) 





/I / /I o\ 

-2.34(48) 


nj 


^ A° 


0.77(10) 


1.57(32) 





2.34(48) 


E+7 


^E*+ 


-0.05(37) 





-0.36(10) 


-0.41(45) 


EO7 


^E*o 


-0.02(17) 


0.01(8) 


-0.32(9) 


-0.31(16) 


E-7 


-» E*- 





0.11(20) 


-0.43(10) 


-0.32(12) 


=07 


_^ w*o 


0.15(24) 





0.23(6) 


0.38(29) 


H"7 







-0.08(12) 


0.26(6) 


0.18(7) 


A7 ■ 


E*o 


0.42(11) 


0.20(6) 





0.62(17) 



D. Individual quark contributions 

To gain a deeper understanding of the dynamics, we 
consider the individual quark sector contributions to the 
transition amplitudes. In our approach, we can easily dial 
individual quark contributions to the QCD sum rules. 
For example, to turn off all u-quark (d-quark) contribu- 
tions, we set the charge factor e„ = (e^ = 0). To turn 
off all s-quark contributions, we set = 0, rus = 0, 
/ = 1, and 0=1. We can extract a number correspond- 
ing to each quark contribution from the slope of Eq. (pD]) 
as a function of 1/M^, while keeping other factors the 
same (/?, Borel window, w). We call this the raw indi- 
vidual quark contributions to the transition amplitudes. 
Table IVl gives the results for the Gi amplitude from such 
a study. We see that the u-quark contribution in the 
transition p"f — > is twice the d-quark contribution as 
expected, and both contributions are positive. It is the 
opposite in the transition —> A*^. The s-quark con- 
tribution is positive in the E transitions, and negative in 
the S and A transitions. In the 'E'^j channel, the 

d-quark and s-quark contributions largely cancel. 

Table|V]gives the results for the G2 amplitude. Similar 
pattern is observed in the proton and neutron channels, 
albeit the signs are opposite for Gi and G2. The s-quark 
contribution is opposite to that in Gi : negative in the E 
channels and positive in the S channels. It donimates in 
the S channels over the u-quark and d-quark contribu- 
tions. Interestingly, s-quark contribution is exactly zero 
in the A channel. 



E. Comparison with other theoretical calculations 

The determination of N to A transition amplitudes 
has been made in a number of other theoretical ap- 
proaches. In the following we briefly discuss each of these 
approaches in order to put our results from QCD sum 
rules in context. The discussion is by no means exhaus- 
tive, but indicative of the breadth of interest in these 
transition amplitudes. First, we present our results in 
Table Ivni after converting from Gi from structure WEi 



and G2 from WEg into the various conventions discussed 
earlier, and compare them with those from a lattice QCD 
calculation [l| , a quark model calculation , and exper- 
iment [l2| . For Ga/1 , we see that our results are slightly 
higher than those from lattice QCD and quark model 
calculations, but the overall pattern (in terms of mag- 
nitude and sign) is consistent among the three very dif- 
ferent calculations. For the ratio Rem, the results are 
very different. We predict negative values for all transi- 
tions at the level of 10% uncertainty. The lattice QCD 
results have too large errors to resolve the sign, although 
a more recent calculation in the — > A"*" channel gives 
a negative value Q . The quark model results have vary- 
ing signs. It demonstrates the difficulty of quantifying 
the small deformation from spherical symmetry in these 
transitions. 

Since pry — )■ A+ is the channel most widely studied 
and the only one measured in experiments, we give a 
more detailed comparison including more theoretical de- 
terminations. Fig. [7] summarizes the calculations of the 
magnetic dipole transition amplitude Gmi hi units of 
nucleon magneton. The calculations are grouped into 
six catagories: lattice QCD (Latt.), QCD sum rules re- 
sults from this work (QCDSR), light-cone QCD sum rules 
(LCSR), hedgehog models including Skyrme and hybrid 
models (Skyrme), quark model calculations (Q.M), and 
bag models (Bag). The experimental result from the 
PDG is also displayed (Expt.). Fig. [5] summarizes the 
calculations of the ratio Rem in the same channel. 

There are only a few lattice QCD calculations. The 
original one by Leinweber et. al. [l| obtained good sig- 
nals for Gmi, but barely a signal for Rem (+3 ± 8 %), 
and no signal for Gc2- Shown here are the more recent 
calculations by Alexandrou et. al. which obtained rea- 
sonable signals for all three form factors. For Gmi, from 
top down, in units of nucleon magneton {^n), the results 
are 2.5±0.4 ,lj with a quenched Wilson action; 1.68±0.06 
with a hybrid action, 1.48 ±0.04 with a quenched Wilson 
action, and 1.07 ± 0.04 with a 2-flavor Wilson action [3|. 
Strictly speaking, we should compare results at the chiral 
limit and — 0. The lattice results quoted are at pion 
mass of around 400 MeV and non-zero momentum trans- 
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TABLE VII: A comparison of N to A transition amplitudes from various approaches: this work (QCDSR) in three different 
conventions; lattice QCD l]; quark model [1^, and experiment 12]. 



Transition 












Lattice QCD 


Quark Model 


Exp eriment 








J Ml 


^1/2 




Gm\ 


Rem 


Gmi 


Rem 


(j-Ml riEM 






/o 


vjre V 


vj»e V 


vjrC V 




% 




% 


/iB /O 


— 

P7 — > A 




-i.DD(i 1 ) 


0.4UO 


-O.iy 


n OCX 
-0.36 


2.46(43) 


3(8) 


2.15 


-0.009 


6.U2 -2.o(oj 


nj — » a'' 


-4 '^Mi'W 


-1.66 (17) 


-0.405 


0.19 


0.36 


-2.46(43) 


3(8) 


-2.15 


-0.009 




E+7 E*+ 


4.18 (42) 


-2.85(29) 


0.240 


-0.11 


-0.22 


2.61(35) 


5(6) 


2.61 


-0.210 




EO7 E*° 


1.79 (18) 


-2.30 (23) 


0.104 


-0.05 


-0.09 


1.07(13) 


4(6) 


1.1 


0.192 




E~7 ^ E*- 


-0.53 (5) 


-7.99 (8) 


-0.031 


0.01 


0.03 


-0.47(9) 


8(4) 


-0.4 


0.985 




=07^5*° 


-1.69(17) 


-1.59(16) 


-0.094 


0.04 


0.08 


-2.77(31) 2.4(2.7) 


-2.86 


0.031 






0.19 (2) 


-12.42 (13) 


0.010 


0.00 


-0.01 


0.47(8) 


7.4(3.0) 


0.44 


-0.259 




A7 ^ E*o 


3.34(34) 


-4.62(48) 


0.314 


-0.135 


-0.285 
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Expt. Latt. QCDSR LCSR Skynne Q.M, Bag 

FIG. 7: A comparison of the magnetic dipole transition am- 
plitude Gmi in units of nucleon magnetons in the p7 A+ 
transition from different approaches. See text for discussions. 
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FIG. 8: A comparison of the Rem ratio in the p7 — > A 
transition. See text for discussions. 



fer. For the ratio Rem, the results are —4. 1 ± 1.4% 
at = 0.13 GeV^ extrapolated to the chiral limit in 
quenched QCD; and -0.77 ± 0.88% H at = o.04 
GeV^ and = 360 MeV in full QCD with a hybrid 
action. The lattice signals for the transition form factors 
are by now well establsihed. With increasing computing 
power and better lattice technology, lattice calculations 
offer the best promise for precision determination of the 
transitions from QCD. 

The two numbers for QCDSR are from different sum 
rules: the top one correponds to Gi from WEi and G2 
from WEg, the bottom one correponds to Gi from WO2 
and G2 from WEg. They are converted to Gmi and Rem 
for comparison. 

The ligh-cone QCD results come from Ref. \6] which 
involve photon distribution amplitudes (light-cone wave- 
functions) up to twist 4. 

The hedgehog models from top down include the holo- 
graphic QCD calculation by Grigoryan, Lee and Yee |^ 
which can be considered as a 5D version of the Skyrme 
model; the SU(3) Skyrme model calculation by Chem- 
tob [3^; the SU(2) Skyrme model calculations by Kunz 
and Mulders [3a|; and the SU(2) Skyrme model calcu- 
lation by Adkins, Nappi and Witten [33|- The semiclas- 
sical hedgehog calculation of Cohen and Broniowski jsl] 
gives a negative value (-2.8) for Ml moment so it is not 
shown in Fig. [71 Skyrme models produce reasonable re- 
sults for the Ml amplitude, but vanishing ones for eletric 
quadrupole E2 transition matrix elements so they are left 
out in Fig. [S] 

The quark model results include (from top down): a 
calculation by Franklin [3^ (It has been correlated with 
Ga and provides arguably the best constraint on Ga/i in 
quark- model calculations); a Bethe-Salpeter determina- 
tion from Mitra and Mittal jl^ ; a calculation by Guiasu 



and Koniuk 41| in which mesonic dressings of the nu- 
cleon are explicitly included; a calculation by Capstick 
(4^ in which configuration mixing in the baryon SU(6) 
wave functions is accounted for; and a calculation by 
Darewych et. al. p3| based on the simple quark model. 
Quark model calculations generally give small values for 
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Rem- 

The bag models include a chiral bag calculation (top) 
by Kalbermann and Eisenberg [3], and an old MIT bag 
calculation (bottom) by Donoghue et. al. [4^ . 

VI. CONCLUSION 

We have carried out a comprehensive study of the 
7V7 A transition amplitudes using the method of QCD 
sum rules. We derived a new, complete set of QCD sum 
rules using generalized interpolating fields and examined 
them by a Monte-Carlo analysis. Here is a summary of 
our findings. 

We proposed a new way of extracting the transition 
amplitudes from the slope of straight lines as a function 
of 1/M^ in conjunction with the corresponding mass sum 
rules, as defined in Eq. ([SO)) . We find that this method 
is more robust than from looking for 'flatness' as a func- 
tion of Borel mass. The linearity displayed from the OPE 
side matches well with the phenomenological side in most 
cases. The method also demonstrates clearly the role 
of the non-diagonal transition terms in the intermediate 
states caused by the external field: wherever such tran- 
sitions are large, the corresponding sum rules perform 
poorly. 

Of the 96 sum rules we derived (from 12 independent 
structures, each for 8 transitions), we find that the sum 
rules from the WEi and WO2 structures are the most re- 
liable for the transition amplitude Gi , based on OPE con- 
vergence and ground-state pole dominance, and small- 
ness of the non-diagonal transitions. The QCD sum rules 
from these structures are in Eq. (|40|) and Eq. (|53|) : their 
predictions are found in Table |lT] and Table IIIII Our 
attempt to extract the G2 amplitude was less success- 
ful. The only sum rules that can give stable results are 
in Eq. (P0|) from structure WEg. Their predictions are 
found in Table IIVI Our final results in the various con- 
ventions are found in Table IVIII 



Our Monte-Carlo analysis revealed that there is an un- 
certainty on the level of 10% in the transition ampli- 
tudes if we assign 10% uncertainty in the QCD input 
parameters. It goes up to about 30% if we adopt the 
conservative assignments that have a wide range of un- 
certainties in Ref. [s^]. The Monte-Carlo analysis also 
revealed some correlations between the input and output 
parameters. The most sensitive is the vacuum suscepti- 
bility X- So a better determination of this parameter can 
help improve the accuracy on the transition amplitudes 
and other quantities computed from the same method. 
We also isolated the individual quark contributions to 
the transition amplitudes. These contributions provide 
insight into the effects of SU(3)-flavor symmetry break- 
ings in the strange quark, and environment sensitivity of 
quarks in different baryons. 

We compared our results with a variety of theoretical 
calculations, and with experiment in the proton chan- 
nel. Our result for transition amplitude Gmi is larger 
than the experiment and other calculations, while our 
result for the ratio Rem is consistent with experiment. 
In general, we find that Gi amplitudes are more stable 
than G2 amplitudes in this approach. Our results for the 
A7 — > transition are new. 

Overall, the results for the N to A transition ampli- 
tudes in the QCD sum rule approach are not as robust 
as those for the baryon magnetic moments [l3, 113, HI] . 
It is in large part due to the intrinsic un-equal mass dou- 
ble pole in the phenomenological representation of the 
spectral functions. Nonetheless, the calculations offer 
a physically-transparent, QCD-based perspective on the 
transitions in terms of quarks, gluons and vacuum con- 
densates. 
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